Quantum Capillary Waves at the Superfluid— Mott Insulator Interface 
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We discuss quantum fluctuations of tlie interface between a superfluid and a Mott-insulating 
state of ultracold atoms in a trap. The fluctuations of the boundary are due to a new type of 
surface modes, whose spectrum is similar — but not identical — to classical capillary waves. The 
corresponding quantum capillary length sets the scale for the penetration of the superfluid into the 
Mott-insulating regime by the proximity effect and may be on the order of several lattice spacings. 
It determines the typical magnitude of the interface width due to quantum fluctuations, which may 
be inferred from single site imaging of ultracold atoms in an optical lattice. 
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The study of fluctuating interfaces is one of the central 
topics in statistical physics |T], with implications for phe- 
nomena like wetting or capillary forces [2] or the physics 
of biological membranes |3] . The origin of interface fluc- 
tuations in these cases is purely thermal. In recent years, 
a lot of interest has focussed on phase transitions that are 
driven by quantum rather than thermal fluctuations [Ij. 
Somewhat surprisingly, the issue of interface fluctuations 
at the boundary between ground states with different or- 
der has not received much attention so far. In our present 
work, we study the interface between a superfluid (SF) 
and a Mott-insulator (MI) realized for ultracold bosons 
in an optical lattice |5] as an elementary example of a 
quantum interface problem. Since the MI-SF transition 
is of second order, the coexistence of ground states with 
different order in this case is due to the presence of a 
trapping potential, which gives rise to a wedding cake 
structure of successive superfluid and Mott-insulating do- 
mains 12]. This nontrivial spatial structure has been ob- 
served in a direct manner recently by a quantum gas 
microscope 13 [5], which provides single site resolution 
of individual atoms in an optical lattice. Within the 
standard local density approximation (LDA), the SF- 
MI interface is sharp. Specifically, for a 2D gas in an 
isotropic trap, it is a perfect circular line. Its position is 
determined by the condition that the local value of the 
chemical potential is equal to the critical value for the 
generic, density driven SF-MI transition of the homoge- 
neous system [51 (see Fig. 1). As we will show below, 
a calculation which incorporates fluctuations around a 
spatially varying smooth background proflle of the su- 
perfluid order parameter near the SF-MI boundary gives 
rise to fluctuations of this interface. They lead to a quan- 
tum uncertainty in its position which can be described 
in terms of an effective capillary length Xg. The spec- 
trum Lu{k) of the elementary excitations, which are lo- 
calized near the interface, crosses over from a gravity 
wave like form uj(k) = ^/g^gk at small wave numbers 
kXg ^ 1 to a free particle like dispersion a;(fc) ~ fc^ 
at kXg S> 1. This is reminiscent of classical capillary 
waves, where the role of gravity is played by the exter- 



nal trap potential. The fc^ behavior at short wavelengths 
is due to the fact that the SF order parameter vanishes 
exponentially as one moves into the MI region. Interac- 
tions between the mobile particles thus become negligi- 
ble. The resulting free particle dispersion is quite differ- 
ent from the k^^"^ behavior found for standard capillary 
waves, which is due to a non-zero surface tension. In 
the SF-MI case, the latter is zero, however, because the 
transition is continuous [HIS]. An important feature of 
this spectrum is that the amplitude of quantum fluctu- 
ations of the position of the SF-MI boundary diverges 
as the density gradient goes to zero. On the qualitative 
level mentioned above, the waves are similar to crystal- 
lization waves at the rough superfluid-crystal boundary 
of ^He [TU]. However, there are important differences be- 
tween these cases. The superfluid-crystal transition in 
''He is of flrst order, with a jump of the density at the 
boundary. As a result the crystallization wave spectrum 
has a form u) ^ fc^/^, and quantum fluctuations of the 
boundary position do not diverge. Moreover, recent ex- 
periments indicate that the superfluid-crystal boundary 
of ^He is quantum smooth and that crystallization waves 
occur only at flnite temperature 



For a quantitative description of interface fluctuations, 
we use an effective action approach, with (j){x, r) the com- 
plex scalar order parameter of the superfluid ^ . At the 
mean field level, the SF-MI transition appears when the 
dimensionless coefficient r of the quadratic contribution 
r|0p to the effective Lagrange density vanishes. In the 
presence of the trap, r{x) is spatially dependent, van- 
ishing at a sharp boundary within LDA. Choosing a co- 
ordinate system where this boundary coincides with the 
y-axis, we have r — bx + . . . locally. The coefficient b 
is determined by the associated gradient of the chemical 
potential fj,. Since all other coefficients of the effective 
action S[(j)] for the order parameter are finite near this 
boundary, the relevant model to describe interface flue- 
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FIG. 1: Qualitative zero temperature phase diagram of the 
homogeneous Bose- Hubbard model. In a trap, the local value 
of the chemical potential varies from a maximum in the center 
to the value where the density vanishes (vertical line). A SF- 
MI interface appears in the vicinity of a generic transition 
point (marked with black dots). Specifically, we will discuss 
the transition on the lower side of the Mott lobe marked by 
the heavy dot. 



tuations near the SF-MI transition is given by 

J d^xdf{e,\y^{x,T)f+bx\cP{x,T)f 

+u\(I){x,t)\'^ + d(j)*{x,T)df(P{x,T)} . (1) 

Quite generally, the coefficients £^Q,b,u and d are phe- 
nomenological parameters. Within a Bose-Hubbard 
model description of the SF-MI transition, they can be 
calculated directly from the hopping and interaction en- 
ergy parameters J and U of the microscopic Hamiltonian, 
using the dimensionless function |12l ITB] 
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Here, g = 1,2,... is the integer density corresponding to 
a specific Mott lobe, while /i = ^'/AJ and U = U/AJ. 
One generally has b = —Ad/^xodxfJ' and d — —dpXo- The 
general expression for u is quite unwieldy, but reduces 
to M « in the regime U ^ J. Note that for the 
lower boundary of the one-atom Mott lobe, (see Fig. 1) 
dfiXo < . Moreover, f = Jr is the dimensionless time, 
while the scale for u and = ^ is simply the lattice 
spacing. For J/U close to the critical value near the tip 
of the lobe, b and d vanish so that higher order terms 
in X must be taken into account for b while the d term 
must be supplemented with the second time derivative 
term following from the gradient expansion of the ac- 
tion. Quite generally, the characteristic length de- 
fines a bare length scale such that the bulk correlation 
length is ^ = ^o/y/r in a homogeneous system, diverg- 
ing as (/i — fj,c)~^^'^ at the transition. Note that we are 



considering the problem in two dimensions, which is the 
upper critical dimension since the dynamical exponent is 
z = 2 for the generic SF-MI transition [11 [9] . As a re- 
sult, mean field theory correctly describes the divergence 
of the correlation length up to logarithmic corrections. 
In the presence of an external trap potential, the co- 
efficient r vanishes linearly in the vicinity of the sharp 
SF-MI interface that results within LDA. The character- 
istic length scale over which this sharp profile will be 
smeared out by fiuctuations is determined by the condi- 
tion S,o/ \/r{Xg) = Xg. It identifies Xg with the scale at 
which the local correlation length reaches Xg itself |14j . 
This results in a broadening of the interface over a scale 
Xg = i£.o/by^'^, a result that is borne out in detail by our 
calculation below. 

The equilibrium order parameter profile can be ob- 
tained by solving the Euler-Lagrange equation corre- 
sponding to the action which reads 

^"-z0-|0|20 = O , (3) 

where we have defined z — x/Xg and (j) — (^/(jig with 
(jig = \JbXgj2u. This equation does not have a closed- 
form solution, but one readily finds the asymptotic be- 
havior, i.e., ^ %/— z for — z ^ 1. For z ^ 1, the nonlin- 
ear term becomes negligible and Eq. ([3| becomes the Airy 
differential equation which is solved by = Ai(2;). For ar- 
bitrary z, the solution can be obtained numerically. The 
resulting order parameter profile ni^\z) — is shown 
in Fig. [2] and is formally identical to that obtained at 
a SF- vacuum boundary |15ffr7] . Note however that in 
the present case the characteristic length Xg is different. 
In particular, it depends explicitly on the interactions 
through the function xo- Unlike the LDA, the mean field 
solution predicts a smooth transition between insulating 
and superfiuid regions over the length Xg . Since all coef- 
ficients in Eq.(|3| are real, the mean field solution (j)o can 
be chosen real without loss of generality. 

In the following, we consider quantum fiuctuations 
around this mean field order parameter profile. To this 
end, we expand the action ([T]) to second order in devia- 
tions from the mean field solution When 0o is small, 
it is natural to consider fiuctuations of the real and imag- 
inary part of (j), i.e., (f) — (pa + ip + itp. Conversely, where 
00 takes appreciable values it is more natural to take 
into account the U{1) symmetry of the action and con- 
sider fiuctuations of ris and 9, where (j) = ^Jn^e^^ . In both 
cases, fiuctuations take on the form of plane waves par- 
allel to the interface while in the direction perpendicular 
to the interface they form a set of modes which must be 
determined from a solution of the Euler-Lagrange equa- 
tions. The equations for the dimensionless versions of 
(y5 and il) (for ease of notation, we omit the tilde in the 
following) read 

if" - (z + 3n(") + fc2)(p + iuj^ = 

(4) 

- (z + nf^ + h^)^ - = , 
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with ns*^-* = (/)o- In the density-phase representation it 
is convenient to define T = ipod and Af = 6ns/(f>Q which 
obey the coupled equations 



T" - if + k^)r - iioU ^ 



(5) 



where / = d'^.n^^ /2n^^ - {d.nf'^ /2nf^f . Asymptoti- 
cally, / z for z > 1 and / ^ -l/4z^ for -z > 1. The 
function /(z) is shown in Fig. [2j 




FIG. 2: The mean field superfluid density ns"' and the func- 
tion / appearing in the difi'erential equations for T and N as 
functions of in dimensionless units. 



In both representations, we have two coupled sec- 
ond order differential equations with non-constant coef- 
ficients. Analytical solutions may thus be expected only 
in limiting cases. There are two opposite regimes where 
the equations can be solved analytically: firstly, the limit 
where the wavelength of the interface fluctuations is small 
compared to Ag, and secondly, the opposite case where 
the wavelength is much larger than the size of the fluc- 
tuation region. 

The first case turns out to correspond to the limit of 
low superfiuid densities where one can neglect the 
term in the action ([T|. Fluctuations then have the same 
general structure as the mean field solution, i.e.. 



(z) = Ai[z - (cj - fc2)]e'(fc2'-"*) 



(6) 



This fixes the frequency oj up to an additive constant 
uj — k"^ which must be determined by comparison with a 
solution to the complete set of differential equations. 

The second case, the long wavelength limit, is just the 
Thomas-Fermi approximation which is familiar from the 
calculation of hydrodynamic modes in Bose-Einstein con- 
densates [18_. In this limit, one obtains (with 0{y,z) ~ 



and the same equation for the superfluid density fluctu- 
ations. This equation has a complete orthonormal set of 
solutions 



(8) 



9nkiz) = V2ke^' Ki-2kz) 



where L„ are the Laguerre polynomials, n — 0,1,.... 
These modes, which are defined for z < only (the 
Thomas-Fermi approximation to the mean field profile 
vanishes identically for z > 0), correspond to the disper- 
sion law (substituting all constants to obtain a dimen- 
sionful quantity) 



2^g&(2n 
d2 



1) 



-k 



,^2{2n + l)kXg , (9) 



i.e., the solutions can be grouped into "branches" char- 
acterized by the integer n, each branch having a gravity 
wave like \/k dispersion. To make this analogy more ex- 
plicit, one may write the lowest branch as a; = \/gcek, 
with (7eff = 2^g6/d2. An equivalent dispersion relation 
has been derived in the context of the boundary of a 
dilute Bose-Einstein condensate in [19J. 
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ze" + 9' ^{k^x + Lu'^ /2)§ = 



(7) 



FIG. 3: The two lowest branches of the dispersion relation 
uin{k) as obtained from a numerical calculation (solid lines). 
The asymptotic behavior for small and large k is indicated 
with dotted lines. 

From these limiting cases, we see that the dispersion 
relation w{k) starts as w ~ y/k at small k and then grad- 
ually crosses over to a, uj ^ k"^ + const behavior as k 
increases. To see how this crossover happens and to de- 
termine the additive constant for each branch in the large 
k regime, we have numerically solved the coupled sets 
of differential equations Q and ([5| using a matrix Nu- 
merov method |5D]. The result is shown in Fig. [s] The 
additive constants for the lowest two branches are de- 
termined by a least square fit to the upper part of the 
dispersion curves as 0.97(2) and 2.27(9). The results are 
in agreement with a prior numerical solution of the same 
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equations in the context of the boundary of a dilute BEC, 
where the lowest band dispersion was derived using a fi- 
nite difference method [21]. The crossover in the lowest 
band indeed happens at kXg ^ 1 so that Ag plays a role 
analogous to that of the capillary length in the physics 
of surface waves on deep water. Note that this is com- 
pletely different from the results of a previous analysis of 
boundary fluctuations of the SF-MI interface by Mariani 
and Stern [55], who found a large k scaling ij{k) ^ k^^'^ 
based on a model with a phenomenological nonzero sur- 
face tension. 

To estimate the spatial extension of the fluctuation re- 
gion, we introduce the phenomenological height variable 
h{y,t) = \g5n{y,Q,t)/dznf\ The mean square fluctua- 
tions {h"^) would diverge if there was only the k^ part of 
the dispersion, but are made flnite due to the crossover 
to gravity waves w ~ \/k precisely as in the case of clas- 
sical capillary waves. In the long wavelength regime, the 
action can be parametrized using the modes Q and then 
reads 

~ 2u ^ \Snk I \ 2 I \Snk 

k,7n,n ^ ^ ^ / \ ■ / 

(10) 

where a;„ = 2TTn/(3u!g are the bosonic Matsubara frequen- 
cies, rendered dimensionless. Note that the combination 
tpdju of the dimensionful Landau-Ginzburg coefficients 
is dimensionless. This representation allows to calculate 
the variances of and 5n. A lower bound to the total 
fluctuations is then given by the contribution from the 
lowest branch and wavelengths smaller than A^: 

a ('.')< = |^j^^^. (") 

i.e., the fluctuations are infrared convergent so that the 
fluctuation region remains conflned to the mean field 
transition region and the interface is quantum smooth. 

To estimate whether the predicted zero point fluc- 
tuations of the SF-MI interface may be observed with 
current experiments, we take the typical case of *^Rb 
atoms in an optical lattice of period I = 532 nm, a lat- 
tice height of Vb = IQAEr which is close to the transi- 
tion point for the MI-SF transition at unit filling g = 1 
{Er = h? /2m\\a,t is the recoil energy) and a central 
chemical potential /x(0) — 1.1/J.c. With these param- 
eters, the effective capillary length Ag is equal to one 
lattice spacing for an isotropic harmonic trap with fre- 
quency 2tt X 16.3 Hz. Moreover, the characteristic fre- 
quency ujg ~ \J gcfi/Xg where the spectrum crosses over 
from a gravity-like form ^ to the free particle ~ k^ 
regime is about 27r x 9Hz. To see quantum fluctuations 
of the interface requires the temperature to be smaller 
than Swg/fce = 0.43 nK. This is quite challenging to 
reach but appears feasible with novel cooling techniques 
like spin gradient demagnetization, where temperatures 



around 0.35 nK have recently been achieved in a similar 
setup |23]. Note that the regime hiOg > k-eT that is re- 
quired to see quantum fluctuations of the interface is op- 
posite to the standard semiclassical limit hw ^ k-oT that 
is usually considered in the thermodynamics of trapped 
BECs [23]. Moreover, for the parameters above, the char- 
acteristic density at the transition is 0^/^ k, 0.125 so that 

[multiplying by nf'' (0)] the density of mobile holes in the 
transition region is about 0.034 per lattice site. Using 
a quantum gas microscope [3 [S] , the smooth non-LDA 
mean fleld proflle can be measured by averaging over a 
sufficiently large number of images. To probe the dis- 
persion relation, one needs to selectively excite individ- 
ual modes which can equally be achieved thanks to the 
single-site addressability of quantum gas microscopes: by 
modulating, e.g., the lattice depth on the single-site level, 
one can achieve values of 2t: /k ranging from 47ri? (corre- 
sponding to the quadrupole mode, with R the radius of 
the LDA transition circle) down to 21. Since Ag can take 
any value from practically zero to several lattice sites, 
this permits to map out the crossover shown in Fig. [3] 
With the stated parameters, the experimental detection 
of the interface's quantum dynamics is certainly challeng- 
ing, but within reach of current experimental technology. 
Additional flexibility may be gained from the use of non- 
harmonic potentials using phase plates. For example, a 
box-like potential with an added localized strong varia- 
tion permits to have the transition happen at a larger 
radius so that more atoms participate. 

In summary, we have discussed the zero temperature 
quantum fluctuations of the MI-SF interface and found 
that the associated dispersion relation leads to a quantum 
smooth surface. An experimental observation of these 
fluctuations requires very shallow trap potentials, where 
the associated capillary length Ag ~ w~^/^ is at least sev- 
eral lattice spacings. From a more general point of view, 
the interface fluctuations we discussed here are just a par- 
ticular case of the rich physics of interfaces in quantum 
phase transitions. For example, in recent years, a com- 
pletely new type of interfaces has turned into the focus of 
research, in which novel phases appear at the boundary 
between two materials with different ground states. A 
striking example is the appearance of a conducting 2D 
electron gas and even tunable superconductivity at the 
boundary between two insulators [251 126] . 
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